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Abstract 

We prove that smooth solutions of non-ideal (viscous and resistive) incompressible magneto- 
hydrodynamic equations satisfy a stochastic law of flux conservation. This property involves an 
ensemble of surfaces obtained from a given, fixed surface by advecting it backward in time under the 
plasma velocity perturbed with a random white-noise. It is shown that the magnetic flux through 
the fixed surface is equal to the average of the magnetic fluxes through the ensemble of surfaces at 
earlier times. This result is an analogue of the well-known Alfven theorem of ideal MHD and is valid 
for any value of the magnetic Prandtl number. A second stochastic conservation law is shown to hold 
at unit Prandtl number, a random version of the generalized Kelvin theorem derived by Bekenstein- 
Oron for ideal MHD. These stochastic conservation laws are not only shown to be consequences of 
the non-ideal MHD equations, but are proved in fact to be equivalent to those equations. We derive 
similar results for two more refined hydromagnetic models. Hall magnetohydrodynamics and the 
two-fiuid plasma model, still assuming incompressible velocities and isotropic transport coefficients. 
Finally, we use these results to discuss briefly the infinite- Reynolds-number limit of hydromagnetic 
turbulence and to support the conjecture that flux-conservation remains stochastic in that limit. 



1 



I Introduction 



If a plasma is sufficiently collisional, then it can be well-described by fluid-mechanical equations. There 
is a hierarchy of such hydromagnetic models, ranging from standard magnetohydrodynamics (MHD), to 
more refined models such as Hall MHD and the two-fluid model, with separate equations for electron 
and ion fluids [521 13] ■ In all of these fluid models the magnetic fleld lines (or magneto- vortex lines) at 
the limit of infinite conductivity are "frozen-in" to the plasma, as first observed by Alfven 1 . The 
properties of magnetic fields for MHD-type models are closely analogous to the properties of vorticity 
fields for the Navier-Stokes equation of neutral fluids [291 [38]. Thus, the Helmholtz-Kclvin theorem on 
conservation of circulations [T51 [55] has an MHD analogue in the Alfven theorem on flux conservation 
[1] , and the Cauchy formula for vorticity [6] has an analogue in the Lundquist formula for magnetic field 
[27] . These Lagrangian properties of magnetic fields are central to many physical processes in plasmas, 
such as magnetic dynamo and magnetic reconnection. It has been claimed, with some justification, that 
"The most important property of an ideal plasma is flux freezing" [23 (section 3.2). 

The fundamental Lagrangian laws of conservation and line-motion hold exactly only for smooth 
solutions of ideal fluid equations, with zero viscosities and resistivities. Recently, however, it has been 
shown by Constantin and Iyer [9] [191 HH] that the analogous laws of vorticity under ideal Euler evolution 
remain for the viscous Navier-Stokes solution as stochastic laws. To formulate these results, the equations 
for Lagrangian fluid particles advected by the Navier-Stokes velocity field must be perturbed by a 
Gaussian white-noise, with amplitude depending upon the viscosity. A random ensemble of fluid motions 
results, depending upon the realization of the white-noise process. To calculate the fluid circulation on 
a given closed loop in the fluid, the loop is evolved backward in time to obtain a random ensemble 
of loops. The circulation on the given loop is the average over the circulations of the ensemble of 
loops at the earlier time. Not only do Navier-Stokes solutions enjoy this remarkable "stochastic Kelvin 
theorem," or conservation of circulations in the mean, but Constantin and Iyer [HI UHl [22] have shown 
that this property also uniquely characterizes the velocity flelds which satisfy the Navier-Stokes equation. 
Furthermore, vortex-lines at any chosen initial time are "frozen-in" to the stochastic fluid flows and thus 
become themselves stochastic. The resultant (deterministic) vorticity at any point at a later time is the 
average over the random ensemble of vorticity vectors that are advected to that point, stretched and 
tilted, by the stochastic flows. These results provide an intuitive way to understand vortex dynamics, 
both stretching and reconnection, in viscous Navier-Stokes fluids. 
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In this paper, we demonstrate similar stochastic conservation laws and "frozen-in" properties for 
non-ideal (resistive and viscous) plasma fluid models. The proofs of Constantin and Iyer [51 (THl HO] 
exploit the Weber formulation [41] of the incompressible Euler equations, and we shall employ here the 
similar Weber formulations of incompressible hydromagnetic models developed in the work of Ruban 
and Kuznetsov [351 125j . Such Weber formulas are implied by the Hamiltonian structure of the ideal 
fluid models and thus have considerable generality. In the case of standard MHD, the Weber formula of 
fI5\ is mathematically equivalent to the generalized Kelvin theorem of Bekenstein and Oron [2\ , which 
provides a second Lagrangian conservation law in addition to the Alfvcn theorem on flux conservation. 
We shall prove that non-ideal MHD at unit magnetic Prandtl number enjoys analogues of both of these 
results as stochastic conservation laws. For general Prandtl number, we shall show that MHD retains 
at least the stochastic Alfven theorem, or flux-conservation in the average sense. This result may be 
formulated equivalently as a stochastic Lundquist formula, according to which the magnetic field vectors 
are "frozen-in" to the stochastic flows and then ensemble-averaged to yield the resultant magnetic field. 
Similar results shall be established also for two more refined hydromagnetic models, given by the Hall 
MHD equations and the two-fluid model of electrons and ions. 

There are several further directions in which these methods and ideas may be developed. It has 
been shown that the Constantin-Iyer formulation of the incompressible Navier-Stokes equation in fact 
corresponds to a variational principle, a stochastic version of the Hamilton-Maupertuis principle of 
least-action |14j . In that work the stochastic Kelvin theorem was shown to arise from a symmetry of 
the stochastic action under the infinite-dimensional particle-relabelling group. The ideal fluid models 
of hydromagnetics are also Hamiltonian in form, as discussed in [351 1251 [5] and references therein. The 
results of the present paper can be derived from stochastic action principles and, in particular, the 
stochastic version of the Bekenstein-Oron generalized Kelvin theorem for MHD can be shown to arise 
from invariance of the stochastic MHD action under particle-relabelling. These results shall be given 
elsewhere. We shall also extend the main results of this paper in following work 15j to hydromagnetic 
models of compressible fluids with anisotropic transport coefficients, i.e. with differing values of viscosity 
and resistivity in directions longitudinal and transverse to the magnetic field. Such refinements have 
importance in applications, but we confine ourselves here for simplicity to incompressible fluids and 
isotropic coefficients. We shall also make only a few brief remarks, in the conclusion section, about 
turbulent hydromagnetics at very high (kinetic and magnetic) Reynolds numbers [HI [121 [El 11] ■ 
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II Standard Magnetohydrodynamics 



We consider in this section the standard incompressible MHD equations [531 IS] for the velocity field u 
and the magnetic field B in three space-dimensions, written (in cgs units) as: 

ftu + (u-V)u = -Vp+ ^JxB + A u (1) 

^tB = V X (uxB) + A A B (2) 

V • B = V • u = (3) 

Here the mass density p, the kinematic viscosity v — fi/ p and the magnetic diffusivity A = rjc/AiT are all 
assumed to be space-time constants. The electric current J is given by the nonrelativistic approximation 
to Ampere's law as J = ^(VxB). Since the speed of light c cancels in eq.([T|), it is simplifying to set 
c — 1 and it is also convenient to assume that density p = 1. We do both throughout this section. 

As did Constantin-Iyer [9l [191 EO] , we take the flow domain fl to be either the 3-torus T'' or else 
3-dimensional Euclidean space M.^. In the latter case, we require that u,B decay sufficiently at infinity. 
The existence of weak solutions to the above MHD system ([I])-® and the existence, uniqueness, and 
regularity of local-in-time strong solutions are established, for example, in [10l[36]. We shall consider 
here only sufficiently smooth initial data uq,Bo G C'^'"(ri) with fc > 3, where C*'''"(il) for fc > and 
a G (0, 1) is the Banach space of functions that are /c-times differentiable with kth partial-derivatives 
Holder continuous of exponent a. We shall use the fact that there exists for such initial data a unique 
solution of ID)-® with u, B e C([to, i/], ^='"(17), A: > 2, for some T^tf~to> 0. It should be possible 
as in [191 [20] to give a self-contained local-existence result, based upon the fixed-point characterization of 
MHD solutions in the theorems below. However, here we find it simpler to give a more direct argument 
that existing smooth solutions possess the stated stochastic conservation laws. 

II. 1 Unit Magnetic Prandtl Number 

Our first result shall be for the case of unit magnetic Prandtl number, when ly — X. In this case, 
we can prove that there are two stochastic Lagrangian conservation laws for solutions of ([l])-©, one 
corresponding to the Alfven theorem 1 and another corresponding to a generalized Kelvin theorem 
[25l [2]. As in [21 [TH [20j, we shall show that the stochastic conservation laws furthermore uniquely 
characterize the MHD solutions. A precise statement of our result is as follows: 
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Proposition II. 1. Divergence-free fields u,'B E C{[tQ,tf],C'^''^{fl)) satisfy the non-ideal, incompressible 
MHD equations with initial data Uo,Bo E C'^'"(ri) for k > 3 and v = \, iff for all closed, 

rectifiable loops C and for all t E [to, tf] (with A — curl^^B) 



w A{-x,t)-dx ~E (f) Ao(a)'(ia 

Jc Jaict) 

(f) u(x, t)-dx = E ® [uo(a) + Bo(a)xR,(a, t)]-da 



(4) 



(5) 



i(at) 

Here a(x, t) are "back-to-label maps" for stochastic forward flows x(a, t) solving 

dx{a,t) ^u{Si{a.,t),t)dt + V2vdW{t), t>to, x(a,to) = a, (6) 

R*(a, i) is the Lagrangian-history charge density (charge per unit area) satisfying 

9tR*(a,i) = -J(X(a,t),i)(VaS(a,t))-\ t > io, R*(a,to) = 0, (7) 

and E in denotes average over realizations of the Brownian motion W(f) in the SDE 

Before we present the proof of the above proposition, let us make a few explanatory remarks. The first 
result (21) may be re-expressed in terms of a flux integral through a smooth bounding surface S of the 
closed loop C, with C ~ dS. It takes the form of a stochastic Alfven theorem, expressing conservation 
on average of magnetic flux: 



/" B(x,<)-dS(x) ==E [ Bo(a)-dS(a 

Js "'a(S,t) 



i(S,t) 

This result is, in turn, equivalent to a stochastic Lundquist formula for the local magnetic field, 



B(x,i) =E Bo(a).VaX(a,i)|5 



a(x,t) 



(8) 



(9) 



The second law ([5]) is a stochastic Kelvin theorem, expressing conservation on average of generalized 
circulation, analogous to the deterministic result of It is equivalent to a stochastic Weber formula, 
corresponding to the deterministic formula of [35[ I25j , 



u(x, t) = EI 



^a(x,i) (^uo(a) + Bo(a)xR*(a,i) 



a(x,t) 



(10) 



where P denotes the Leray- Hodge projection onto divergence- free vector fields [7]. To explain the physical 
meaning of this formula, it is useful to transform the final term to Eulerian form: 



u(x,t) =EP V^a(x,i)uo(a(x,t)) H- B(x, <) xR(x, i) 



(11) 
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where B(x, i) = Bo(a)-Vax(a, i)!^^^ and R(x, i) = R*(a, <)-Vax(a, t) 



. Then, employing the 



usual rules of calculus, R(x, t) satisfies the Stratonovich SPDE which follows from ([7]): 



dR(x,t) = Vx (U(x,odi)xR(x,t)) - 3ix,t)dt, t > to, 



R(x,to) =0, 



(12) 



with U(x, t) — J^^ dt' u(x, i')+V2i'W(t) the (Ito and Stratonovich) infinitesimal generator of the stochas- 
tic Lagrangian flow. Consider the flux integral of R.^ through any smooth surface: 




Then, (fT2|) is equivalent to the following equation, valid for all smooth surfaces S: 



dQ{S,t) = -~dt- 



L 



Si{S,t) 



J(x,t)-dS(x), t > to, 



Q{S,to) = 0. 



(14) 



The above equation implies that —Q{S,t) equals the electric charge which flowed across the material 
surface x(S', t) between times to and t. This explains the name "Lagrangian-history charge density" for 
the fleld R*(a, i) used in the above proposition. For any inflnitcsimal vector surface element dS(a.) 
starting at point a, — R,(a, t)'(iS(a) equals the charge crossing the advected surface from to to t. 

Proof of Proposition \II. l[ We first remark that, given the velocity field u e C{[to,tf],C''^°'{n)) for fc > 1, 
there exists a stochastic fiow x(a, i) of C*'''"-diffeomorphisms solving the SDE ([6]), so that the inverse 
map a(x, t) exists and belongs to C{[to,tf],C'''°'{il.)), at least for T sufficiently small. This follows by 
the arguments in [T9j [20] or the general methods in the monograph [24 , Chapter 4. If we assume 
that u,B e C{[to,tf],C'''"{n)) for fc > 3, then VaX,J e C{[to,tf],C''^°'{fl)) for fc > 2, so that R, e 
C^{[to,tf], C'^'"(ri)) for fc > 2. This is sufficient regularity to justify all of our calculations below. It is 
particularly important that R, is bounded variation in time and has no martingale part. 

We begin by showing the "if" direction. Therefore, assuming that divergence-free fields u, B e 
C{[to,tf],C'''"{n)) with fc > 3 solve the fixed-point problem (FPP) specified by the cqs.^-jT]) we must 
prove that they also satisfy the incompressible MHD equations ([I])-©. The argument closely follows 
the proof of Theorem 2.2 in Section 4 of We make the successive definitions 



w(x,t)= Uo(a) + Bo(a)xR,(a,t) 



(15) 



v(x,i) = V^a(x,i)w(x,i) 



(16) 



u(x, t) = Pv(x, t) = v(x, t) - V,^(x, t) 



(17) 
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We see that w,v,u e C{[to,tf],C'''°'{Q)) for k > 2 and that the stochastic Weber formula PU)) is 
restated as u = E(u). We now develop a stochastic evolution equation for each of these three variables. 
Note first that the "back-to-labels map" a satisfies 



(18) 



This is proved in [S], Proposition 4.2. It can also be derived as a special case of the "first Ito formula" 
for a backward flow; see [24], Theorem 4.4.5. Then, as a consequence of the generalized Ito rule, 

dw(x,t) = [-(u-Vj;)w + !^ A w]dt+ |"Bo(a)x9tR,(a,i)j dt - V2i^{dW{t)-V (19) 

L J a(x,t) 

For example, see [9], Corollary 4.3. The term in (fT9|) involving 9tR, can be evaluated using 



'V^a{x,t) Bo(a)xafR»(a,t) = J(x, i) X B(x, t), 

L J a(x,t) 



which follows from ([7]), and from the definition 



(20) 



B(x,t) = Bo(a).VaX(a,t)| 



a(x,t) 



(21) 



We next calculate the differential of v using the Ito product rule. 



dv(x, t) = V:ra(x, i)dw(x, t) + d( V:ra)w(x, t) + d{V^a., w) , 



which, together with HI]) , (HH) , ([201) , gives 

dv(x,<)= -(u-Vj,)v - (Va;u)v + J(x,t)xB(x,t) + A V - \/2^(dW(t)-V^)v. (22) 

The rest of the argument goes exactly as in Section 4 of [9]. As in the proof of Theorem 2.2 of [9], the 
differential of u, as defined in p?]) . can be expressed as 



du(x, t) 



-(u-V:,)u - (V:,u)u + J(x,i)xB(x,i) + vAu dt- V2^{dW{t)-Va:)u 
dip + {{u-'Vx)if -iyAip)dt + V2iy{dW{t)-'Vx)'f 



(23) 



and, using again the generalized Ito rule and definition (|2ip . 



dB(x,i)= -(u-V,OB + (B-V^)u + i^AB dt - V2^{dW{t)-Vx)B 



(24) 



just as in the proof of Proposition 2.7 of [9j. Taking the expectation over the Brownian motion in 



eas. (P5)) -([ ^ yields eqs.ld])-© with A = and kinematic pressure p — ^jup + (p + {u-'Vx)'P — v Api. 
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We finally show the "only if" direction. Therefore, assuming that divergence-free fields u. B e 
C{[to, tf], C'^'"{fl)) with fc > 3 solve the incompressible MHD equations ([I])-© we shall show that they 
also satisfy the fixed-point problem (FPP) specified by the eqs.([4])-([7]). Let us define 



u(x, t) = EI 



V:,a(x,i) (uo(a) +Bo(a)xR*(a,t)) 



a(x,t) 



and 



(25) 



(26) 



B(x,t)=E^Bo(a).VaX(a,0|~(, ,)J , 

where x(a, t) solves (jH]) and R*(a, t) solves ([7]), for the given u, B. It then follows from our previous work 
that u, B e C{[to,tf], C'^'"(ri)) with k >2, are divergence-free, and solve the linear equations 



dtu = -(u-V)u- (Vu)u- Vp + JxB + i/Au, 



(27) 



dtB^-{u-V)B + B-Vu + iyAB, (28) 

with initial conditions u(to) — Uo,B(io) — Bq. At least one solution is the pair (u, B) itself, so that, 
if solutions of the initial-value problem are unique, it must be the case that (u, B) — (u, B). It thus 
suffices to prove that the linear system (P7|) - (l^51) has unique solutions for specified initial data. 

This may be shown by a standard argument based on an energy estimate (e.g. see [30]). For the pair 
z(x) = (u(x),B(x)) define norms 

1/2 



Z 2 



d'x [|u(x)p + |B(x)n 



lz||oo = sup[|u(x)|-H|B(x) 



An easy calculation then shows for any solution z(x, i) — (u(x, t), B(x, t)) of (|27ll .([28 | that 



dt 



= 2/ d^a: [c'jUi(x,i) (S,(x,t)Bj(x,t) - M,(x,t)uj(x,i)) 
Jn 

+djB,{x,t) {u,{x,t)Bj{x,t)~Uj{x,t)B,{x,t))] - 2i/||Vz( 
< 2||Vz(0||oo||z(0||^-2H|Vz(i)||i 



d 
dt 



For some e > 0, choose 7 > sup^g^t^j || Vz(<)||oo + £• Then it follows that 

-2^(*-*«)||z(t)||2] < _2e-2^(*-*«) [e\\z{t)\\l + i.\\Vz{t)\\l] . 
Integration yields the energy inequality 

e-27(*/-*o)|,^(i^)||2 + 2 T'dte-^-^^*-*") [6||z(0||^ + H|Vz(t)||^] < ||zo||i 

Jto 

which implies uniqueness of solutions of (P7)) . (P5)) as a direct consequence. 



(29) 
□ 
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II. 2 General Magnetic Prandtl Number 

An examination of the proof in the previous subsection reveals an interesting fact that the equation 
((27)) for u involves both u and B, but the equation ((28)l for B involves only B itself. This implies that 
the stochastic representation previously employed for both u and B can be exploited for B alone and, 
furthermore, at any magnetic Prandtl number. A precise statement of the result is as follows: 

Proposition II. 2. Divergence-free fields u, B G C{[to,tf],C'^''^{fl)) satisfy the non-ideal, incompressible 
MHD equations ([2P-(0) with initial data Uq, Bq G C'^'"(r2) for k > 3 iff the momentum equation {ip holds 
over that interval and simultaneously the stochastic flux conservation holds 



/ B(x,t).dS(x) =E / Bo(a)-dS(a) 

Js Ja(S,t) 



(30) 

5(S,t) 

for all smooth surfaces S and all times t G [to,tf],, where a(x, i) are "back-to-label maps" for stochastic 
forward flows x(a, t) solving the SDE 

c6c(a,t) = u(x(a,t),t)dt + \/2AdW(t), t > io, x(a,to) = a. (31) 

Proof of Provosition UTTBi The argument is nearly the same as that for the previous proposition. The 
"if" direction is immediate, since the result ([50]) is equivalent to the stochastic Lundquist formula Q 
and the equation ^ for B follows from ^ using the generalized Ito rule, just as before. For the "only 
if" direction, we define 

B(x,t)=E[Bo(a).V,S(a,i)|5(^^,) 

where x(a, t) is the stochastic flow defined by the SDE pTjl for the velocity u that satisfies the MHD 
momentum equation ([T]). It follows that B satisfies the kinematic dynamo equation 

dtB = -(u-V)B + B-Vu + A A B, t > io, B(to) = Bq 

one of whose solutions is B = B. Unicity of this solution again follows from an energy inequality 



e-27(t/-to)||B(t^)||2 + 2 / dte-2^(*-*°) [eWmWl + A||VB(i)||2] < ||Bo||i 



which is derived by a similar calculation as before [30j, with 7 > sup(g[jjj (^^j || Vu(t)||oo + e- □ 
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Ill Other Incompressible Plasma Fluid Models 

We now establish similar stochastic conservation laws for some other non-ideal plasma fluid models, more 
refined than standard MHD. Keeping within the stated limitations of this paper, we consider here only 
the versions of these models assuming incompressible fluid velocities and isotropic transport coefficients. 
We also give only sketches of the proofs of the stated theorems, emphasizing essential differences from 
those given previously, since most of the details are very similar. 

III.l Hall Magnetohydrodynamics 

The equations of incompressible Hall magnetohydrodynamics (HMHD) have the form: 

(9tu+ (u- V)u = -Vp+ --^(VxB)xB + i/V^u (32) 

47rp 



dfB = V X 



u- -^VxB 1 X B 



A A B (33) 

V • B = V • u = (34) 

The magnetic induction equation (|33p contains a "Hall drift term" proportional to a = mc/e, whose 
importance was first emphasized by Lighthill 26J and which was subsequently extensively investigated; 
see [H] and p51 [^. The limit a ^ formally recovers standard MHD. Mathematical properties of 
HMHD solutions (existence, regularity, etc.) are studied in [3T[ [32] . 

Before stating our new theorems, we must review some standard facts about HMHD eqs. ([5^ - (|34p . 
for which, for example, see [35j . If one introduces a vector potential A for the magnetic field in Coloumb 
gauge, V-A — 0, then, along with a corresponding scalar potential $, it satisfies 

atA= (u-— VxB) xB- V$ + AAA. (35) 
V 47r/? J 

HMHD is a Hamiltonian fluid model with two canonical momenta 

Pj = u + a~^A, Pe = — Q!~^A, (36) 
which are both divergence-free, V-po- =0, a — i,e. When v — X, these satisfy the equations 

dtPa = UaX{Vxpa) - V-Ka + I' Apa, a^i,e (37) 
with TTi = p + (l/2)|up + a^^^,TTe = — and with the ion fluid velocity and Ug the electron 



fluid velocity, given by 



Ui = u, Ue^u-- — VxB. (38) 
47rp 
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Note that, if 7^ A, then the equation for pi would contain an additional term (i^— A) Apg. Corresponding 
to the two canonical momenta there are two generalized vorticities Jl^ = Vxp(j, c = i, e, or concretely 

ni=Lj + a^^B, J7e = -a^^B. (39) 

When ly — X, these generalized vorticities satisfy 

dtfla = V x{u„xn„) + iyAn„, a^i,e. (40) 

These equations imply two "frozen-in" fields for ideal HMHD, one for the ion fluid and one for the 
electron fluid, and two Cauchy-type formulas for the two generalized vorticities. There are likewise two 
Kelvin-type theorems and two Weber formulas for the two canonical momenta. 

We now state stochastic analogues of these results for non-ideal HMHD, flrst for unit Prandtl number: 

Proposition III. 1.1. Divergence-free fields G C([io, t/], C'^'"(r2)) satisfy the non- ideal, incompress- 
ible HMHD equations i3^) - [FS\) with initial data Uo,Bo G C'^''"(r2), with fc > 3 and for unit magnetic 
Prandtl number vj\— 1, iff for all closed, rectifiable loops C and for all i G [to,tf] 



f pcr(x, t)-dx = E f Pao{a)-da 

Jc Ja„(C,t) 



cr^i,e (41) 



Here the canonical momenta p^-, a — i,e are given by eg. 1136]} and a^{x^t) are "back-to-label maps" for 
stochastic forward flows Xcr(a, t) solving, 

dx.^{a,t) =u„{Si„{ai,t),t)dt + V2iydW{t), t > to, x^(a,to)=a, (42) 

for a = i,e, with velocities given in eg. i38\) and W(t) a standard Brownian motion. 

Remarks: (i) If the noise terms in eq. ((42|) were chosen to be instead \/2va dW{t), a — i,e with Vi ^ v 



and Ve — A, then one would obtain the correct induction equation ([33]) but the momentum equation 
would differ from ([5^ . containing an additional term a~^{v — A) A A = ^^{^ ~ v)i. 
(a) As we shall see below, it would be enough to assume u G C{[to,tf],C'^'°'{Q)) for k > 2, whereas 
fc > 3 is required for B because of the Hall drift term. 

(Hi) The two stochastic Kelvin theorems in (|41|) are equivalent to stochastic Weber formulas 

p^(x,t) =EP[V^S<,(x,t)p^o(S^(x,<))], cr^i,e. (43) 
There are likewise two stochastic Cauchy formulas for the two generalized vorticities: 

f2^(x,i)=E r2^o(a)-VaX^(a,t)|5_^(^_j) , cr = i,e. (44) 
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Sketch of proof of Proposition III. 1.1. : The proof is very similar to those given in the previous section 
and in [9] . The main step is to derive equations for the stochastic time-differential of the variables 

Pa = P[Va;a^(p^o ° a^r)] = ^xaa{Pcro°a„) - Vj;^ (45) 

for a = i, e, of the form 

dp^(x,t) ^ [-{u^-V^)pa-{Va:U^)p^ + ,yAp^]dt~V2^{dW{t)-V^)Pa 

-V^ dipa + {{u^-V^)lpa-vAipa)dt + V2^{dW{t)-V^)lp„ . (46) 

The calculations are essentially identical to those presented before. There is just one technical issue 
related to regularity of the field Pe, which should belong to C^'"(r2) in order to give classical meaning 
to the Laplacian term v A Pe. However, if u,B e C^'°'{Q), then e C^'"(il),Ue e C^'"(f7), so that 
x,,a, e C^'"{n),Xe,ae e C^'°'{n), and thus V^^a, e C'^'°'{n),V^ae G C^^"(f7). This means that Pe 
defined by (|45p belongs a priori only to C^'°'{fl). However, we may use an integration-by-parts identity 

proved as Lemma 3.1 of [20], in order to rewrite the partial derivative dkPe m terms of VajSg only and 
eliminate second derivatives of Sg. We can thus conclude that Pe € C'^'°'{Vl). 
Ensemble-averaging the equations for tr = «,e yields 

diPa = -(U^-V:r)p^ - (V^U<,)p^ - V:,7f^ + V Ap„, (47) 

with Tfff = dt^^ + (U(j-Vi:)^g. — V /\ Tf^, for cr = i, e, or, in terms of u and B variables, 

5tu = -(u • V)u - (Vu)u - Vp-h -^(VxB)xB + A u, (48) 



aB = V X 



u- — VxB 1 X B 



t/AB. (49) 



A fixed point satisfying (u, B) — (u,B) must therefore also obey the HMHD equations ((32| - (|33p . The 
converse statement is obtained from the unicity of solutions to the initial-value problem for the above 
linear equations — either (|T7)) . cr = z, e or (|48p . (|15)) — which is proved using energy estimates as before. □ 

Just as in the case of standard MHD, we see that the equation (|49p for B does not depend upon u. This 
makes it possible to derive a stochastic conservation law for magnetic-flux at any Prandtl number: 
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Proposition III. 1.2. Divergence- free fields u,B G C{[tQ,tf],C^'°'{il)) satisfy the non-ideal, incom- 
pressible HMHD equations !i3S^) - lS^) with initial data Uo,Bo G C^''^{Vl) for k > 4 iff the momentum, 
equation iS^) holds over that interval and simultaneously the stochastic flux conservation holds 



/ B(x,t)-dS(x) =E [ Bo(a)-dS( 

Js Ja(S,t) 



(50) 



i(S,t) 

for all smooth surfaces S and all times t g [t(),tf],, where a(x, i) are "back-to-label maps" for stochastic 
forward flows x(a, t) solving the SDE 

dX(a,<) Ue(X(a,t),i)di + \/2AdW(^), <> to, x(a,to)=a (51) 

with Ue the electron fluid velocity given by h3^). 

The proof is as in Proposition |lL2l but greater smoothness of B is required to guarantee that B g C^'". 
III.2 Two-Fluid Plasma Model 

The most general hydrodynamical model of a fully ionized plasma consisting of electrons and one species 
of singly-charged ions is the two- fluid model of Braginsky [4] , or the Braginsky equations. See also |231 13] . 
The basic variables of the model are the two fluid velocities Ucr, a — i,e, with a = i for the ion fluid 
and <7 — e for the electron fluid. In the simple form considered here the equations take the form: 

(9t -^u,-V)u, = (e+ iu,xB ) - Vpi + Au, - -(u, - Ue), V-u^ = (52) 
■nii \ c ) n 

(a* Ue-V)Ue = ^ (E-H-UeXB ) - Vpe + J^eAUe-— (Ue-Ui), V-Ug = (53) 

rUe \ C J Te 

Ivr Itt 

- AA = — J = — ne(u, -Ue), V-A = (54) 
c c 

F, = --dtA, B = VxA (55) 
c 

In addition to internal viscosities v^i cr = i,e, there are linear drag terms which represent the exchange 
of momentum between the two fluids by collisions of the constituent particles and which are proportional 
to the collision frequencies a ~ i,e. Conservation of momentum requires me /re — rui/Ti. Note 

that the vector potential A is not an independent variable, but is completely determined from Ue, u^ by 
means of the elliptic equation (I54p . Mathematical properties of solutions of these equations (existence, 
regularity, etc.) are studied in [33j . including even a separate equation for neutral molecules. 
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Neglecting viscosities and drag, the two-fluid model is Hamiltonian with canonical momenta: 

6 6 
Pi = H A, Pe = Ue A, (56) 

niiC nieC 
satisfying V-pcr = for cr = i,e. E.g. see [35]. Using (ucr'V)ucr = V(^|uctP) — UcrX(VxUcr), i^a = 
Pa + -^luo-p, the eas. (|5^ .(|55 1) for Ue, can be rewritten for Pe, pi as 

dtPi = UiX(Vxpi) - VvTi + v^A\ii - — (ui - Ue). (57) 

Ti 

dtPe = UeX(VXPe) - VlTe + I'e A Ug - — (Ug - Ui), (58) 

Define the magnetic diffusivity A = rriiC^ / A-Kne^Ti — rrieC^ /Anne'^Te, so that ^(ug — u^) = ^(ug — U;) = 
A • I A A. Hence, choosing i^g = z/^ = A, ([57]) . ((55|) become 

StPo- = u^x(VxPct) - Vtt^ + A A Pct, (T = i,e. (59) 

The vector potential A can be recovered from p^, Pe by solving the Helmholtz equation 

- AA + K^A = — ne{p, - Pe), V-A = (60) 

c 

with = ATTne^ / iic^ and — m^^ + m^^ and then Ui,Ue obtained from (|56p . For this non-ideal 
version of the two-fluid model there are two stochastic conservation laws corresponding to the two 
canonical momenta: 

Proposition III. 2. Diverg ence- free fields Ue,Ui G C([to, ^/], C'^'"(f^)) satisfy the non-ideal, incompress- 
ible two-fluid equations I5ii|) -i f55)j with initial data Ueo,Uio S C'^'"(ri) for fc > 2 and for unit magnetic 
Prandtl numbers v^jX — = 1, iff for all closed, rectifiable loops C and for all t G [to, tf] 



(p Pct(x, t)-dx = E J) Po-o(a)'rfa 



a = i,e (61) 

Here the canonical momenta p^, a ~ i,e are given by eg. I156\) anda.a{^,t) cire "back-to-label maps" for 
stochastic forward flows Xo.(a, t) solving, 

dXa{Bi,t) = Ua{xa{a,t),t)dt + V2XdW{t), t> to, X„(a,io)=a, (62) 

for cr = i,e, with W(/;) a standard Brownian motion. 

Remarks: (i) li the noise terms in eq. (p^ were chosen to be \J2\c (T^ ait), a — i,e with Ae 7^ \i then 
one would obtain two-fluid model (|52p - (|55p with = Xa and = ma-c^ / Anne'^ for a — i,e. Although 
mathematically well-posed, this system is unphysical since it violates conservation of momentum. 
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(a) The two stochastic Kelvin theorems (j61|) are mathematically equivalent to stochastic Weber formulas: 

Pff(x,t) = EP[V:rao.(x,t)p^o(aff(x, i))] , a = i,e, (63) 
identical in form to (HI) for HMHD. 



Sketch of proof of Provosition [Til. Si The proof is very similar to that of Proposition III. 1.1 and to the 
proofs in Section 4 of [S]. Stochastic Weber variables Pa, (t — i,e oi the same form as (|45l) arc shown to 
obey stochastic PDE's of the same form as pS)) . It is now enough to assume Ug, e ^([to, tf], C*'''"(il)) 
for fc > 2, because the integration-parts-identity [19l[20] can be employed to show that Pe, Pi G C^'"{n)), 
just as for pi in the proof of Proposition III. 1.1. 

□ 

The curl of the two canonical momenta in (|56p give two generalized vorticities, fla = Vxpo- for a — i,e: 

=U3i + — B, = We ^B. (64) 

rriiC rrieC 

If we assume sufficient smoothness (Ug, G C{[to,tf], C'^'°'{il)) with k > 3), then these satisfy equations 

dtn^ = V x{u„xn„) + XAn^, a^i,e, (65) 

and there are two stochastic Cauchy formulas 

r2ff(x,<)=E O^o(a)-Vax^(a, i)|5^(^ J) , a^i,e, (66) 

of the same form as (|33]) for HMHD. Unlike for the previous models, for the two-fluid model there is no 
separate stochastic "frozen-in" property for the magnetic field B alone, at general Prandtl number. 



IV Discussion 

The results of the present paper provide new tools with which to investigate and explain resistive 
phenomena in plasma fluids. The stochastic Lagrangian conservation laws derived here are the analogues 
for non-ideal hydromagnetic systems of flux-freezing for ideal ones. Especially robust is the stochastic 
Alfvcn theorem and stochastic Lundquist formula, which hold in both MHD and HMHD at any magnetic 
Prandtl number. These results have important implications for resistive magnetic reconnection and 
related problems such as magnetic dynamo [23l[3], which will be pursued in detail in future publications. 
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Here we just note that that the stochastic Lundquist formula describes how the resuhant magnetic field 
B(x, t) at a spacetime point (x, t) is obtained by advecting all magnetic field lines as "fi'ozen-in" to the 
stochastic flows, with added white-noise, and then averaging those magnetic field vectors that arrive 
to the given point. The advection by the stochastic fiows produces the nonlinear eff'ects of magnetic 
stretching and tilting, while the average over the Brownian motions represents the resistive "gluing" 
of the magnetic field, reconnecting the field-lines and changing their topology. Such resistive effects 
are recognized as important in the dynamo process by the cycle of "stretch- twist- fold- reconnect" [17) . 
Moffatt has referred to the "oxymoronic role" of resistivity, writing that "the dynamo process may be 
described as a process of 'regenerative decay', or perhaps better 'reinvigorating dissipation'." [28j 

A possible criticism of the physical relevance of our results is that molecular resistivity, represented 
by a Laplacian term in the induction equation, is a poor model of actual dissipative processes in a 
plasma. In contrast to the iconic status of the viscosity term in the Navier-Stokes equation for neutral 
fiuids, there is considerably less universality in the form of the dissipation in plasmas or, indeed, in the 
validity of a hydromagnetic description. The two-fiuid equations of Braginsky |4j (see also ^51 [5])are 
more complicated than those discussed in our section [111.21 For example, viscosity and resistivity in the 
standard Braginsky equations are anisotropic, with magnitudes differing along directions longitudinal 
and transverse to the local magnetic field. Furthermore, microscopic Spitzer resistivity is not the only 
form of magnetic dissipation that may occur in plasmas. A wide variety of processses, both coUisional 
and non-coUisional, have been proposed to lead to "anomalous resistivity" of different forms [34l [40] . 
Furthermore, in a partially ionized plasma the collisions of ions with neutral molecules induces an 
"ambipolar drift" of magnetic field lines with velocity proportional to the Lorenz force [371 [S] and this 
can be the most significant form of magnetic dissipation in some cases, e.g. the interstellar medium. 
Thus, the fluid models that we have considered are not necessarily the most physically realistic. 

There are two responses that we can give to this important set of criticisms. 

First, the results presented in this paper are far from the most general possible. We have chosen to 
restrict discussion here to models with incompressible fiuids and isotropic transport coefficients, since 
these hydromagnetic models are widely employed and the proofs of the main results are simpler for 
them than for more complete models. However, in a following work [15j . we establish similar results 
for much more general plasma fluid models, allowing for compressible fluids, anisotropic pressure and 
transport coefficients, neutral components, etc. Stochastic conservation laws of the sort demonstrated 
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here are quite general and should hold for a very large class of non- ideal plasma fluid models, when the 
ideal version of the model possesses a corresponding "frozen-in" field. It may even be possible to prove 
similar stochastic laws for kinetic models of plasmas with collisions described by Boltzmann kernels |16| 
or Fokker-Planck operators [23] (Section 8.3), since the ideal, coUisionless Vlasov dynamics possesses 
analogues of the frozen-in invariants |43j . 

Second, the precise form of the dissipation in hydromagnetic models may not matter, as long as its 
effects are confined to sufficiently small length-scales. There is then a large "effective Reynolds number" 
(both magnetic and kinetic) and the plasma fluid becomes turbulent. We have previously argued [12] 
that the laws of flux conservation and magnetic line-motion in hydromagnetic turbulence are intrinsically 
stochastic in the limit of infinite Reynolds number. Formally, the random white-noise disappears in the 
equations for stochastic Lagrangian particles 

dx = u''{x{t),t)dt + V2XdW{t), t>to, x(to) = xq (67) 

as i^, A — > (cf. also eq.®). However, the randomness need not vanish if the advecting velocity 
u"^ solving ([T]) approaches a rough or singular velocity u in this limit, as expected for a Kolmogorov- 
type cascade range. As a consequence of "explosive" separation of particles in Richardson two-particle 
turbulent diffusion, a pair of solutions of (j67|) with the same initial condition xq may separate at time t 
to a mean-square distance ~ in the limit u, X 0. These statements have been proved rigorously to 
hold in the Kazantsev-Kraichnan kinematic dynamo model [HJ [52] . It has furthermore been proved that 
the Lagrangian trajectories in the Kazantsev-Kraichnan model remain stochastic as i/, A — > 0, a result 
that has been termed "spontaneous stochasticity." The limiting probability distributions of trajectories 
are known to be very robust and universal for the case of an incompressible fluid velocity, with the same 
result being obtained for limits of a wide class of regularizations. See |111 112] for references and more 
detailed discussion. The rigorous results for the Kazantsev-Kraichnan dynamo model and the new results 
in the present work give further plausibility to the ideas that the precise form of dissipation docs not 
matter in nonlinear hydromagnetic turbulence and that flux conservation and "frozen-in" line-motion 
will remain as stochastic laws in the limit of very large Reynolds numbers. 
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